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ABSTRACT 



Tensor-scalar theory of gravity allows the generation of gravitational waves from astrophysical 
sources, like Supernovas, even in the spherical case. That motivated us to study the collapse of 
a degenerate stellar core, within tensor-scalar gravity, leading to the formation of a neutron star 
through a bounce and the formation of a shock. We discuss in this paper the effects of the scalar 
field on the evolution of the system, as well as the appearance of strong non-perturbative effects 
of this scalar field (the so-called "spontaneous sealarization" ). As a main result, we describe the 
resulting gravitational monopolar radiation (form and amplitude) and discuss the possibility of 
its detection by the gravitational detectors currently under constnietion, taking into account 
the existing constraints on the scalar field. From the numerical point of view it is worthy to 
point out that we have developed a combined code which uses pseudo-spectral methods, for the 
evolution of the scalar field, and High Resolution Shock-Capturing schemes, for the evolution of 
the hydro dynamical system. Although this code has been used to integrate the field equations of 
that theory of gravity, in the spherically symmetric case, a by-product of the present work is to 
gain experience for an ulterior extension to multidimensional problems in Numerical Relativity 
of such numerical strategy. 

Subject headings: Gravitation - Gravitational waves - Methods: numerical - 
Supernova: general - Shock waves 
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Introduction 



Tensor-scalar theories of gravitation (Damour & Esposito-Farese (1992)) appear to be interesting 
alternatives to Einstein's General Relativity, which naturally arise from the low-energy limit of superstring 
theories ( Damour fc Polyakov (1994)[ ). These theories also allow for a "graceful exit" from inflation 
( iGarcfa-Bellido fc Quiros (1990)| ). One of their main observational predictions is the emission of monopolar 
gravitational radiation, from quasi-spherical astrophysical objects. Previous works on collapses in 
tensor-scalar theories dealt with the gravitational collapse in tensor-scalar theories started with studying 



the Oppenheimer-Snyder (dust) collapse (Oppenheimer & Snyder(1939)) in the frame of the so-called 



"Jordan-Fierz-Brans-Dicke" theory, in which the coupling function (^) depends only on one parameter 



(Shibata et al. (1994) and Scheel et al. (1995)). The next work by Harada et al. (1997) considered a 
two-parameter coupling function (as in this work), but neglecting any non-perturbative effect from the 
scalar field. However, these non-perturbative effects can be very important in neutron stars, as have shown 



Damour fc Esposito-Farese (1993) , making an analogy between this "spontaneous scalarization" of neutron 



stars and the spontaneous magnetization of ferromagnets ( Damour fc Esposito-Farese (1996) ). Novak 
(1998a)| (paper I hereafter) dealt with the collapse of a neutron star to a black hole, including pressure 



and non-perturbative effects. We here study the emission of monopolar gravitational waves coming from 
the formation of a neutron star. The problem is quite different since, according to Hawking's "no-hair" 



theorem (Hawking (1972)), a black hole cannot possess any scalar charge, neutron stars appear to be the 
astrophysical object with the strongest possible scalar charge. Their formation (during a Supernova type 
lb or type II event) is therefore a powerful and interesting source of gravitational radiation. The purpose 
of this paper is to present results from the numerical study of the collapse of a stellar core which forms a 
neutron star, within the framework of tensor-scalar theories and, namely, to explicit the resulting monopolar 
gravitational radiation. 

Following paper I, we will consider only one scalar field ip and we will use two conformally related 
metrics g^j, = a?{ip)g*^^, called the Jordan-Fierz (or physical) metric {g^v) and the Einstein metric (g^j,). 
All the quantities with a tilde are related to the Jordan-Fierz metric, to which the matter is coupled. That 
means that all non-gravitational experiments measure this metric, although the field equations of the theory 
are better formulated in the Einstein one. The tensor-scalar field equations are 



R 



1 



SttG, 



(1) 



-4- 



for the tensor part, and 

5rv;V> = -4^G,a(v5)r, (2) 

for the scalar field. 

The Einstein-frame stress-energy tensor T^'^ is related to the physical one by the arbitrary coupling 
function 0(1^9) 

Tf, = a\ip)f{;. (3) 
The function a(ip) is the logarithmic derivative of a(ip) and represents the field-dependent coupling strength 



between matter and the scalar field. We will assume it is well represented by a linear function (see Novak 
(1998b)| for more details): 

a{ip) = ao + Po X {ip - Lpo). (4) 

The arbitrary function a{ip) is determined by the two parameters (ao,/3o) ivo^ the scalar field at spatial 
infinity, is redundant with ao). We will also consider the collapse to be spherically symmetric, since we are 
concerned with monopolar radiation, which should be dominant over higher terms in such events. 

The paper is organized as follows: Section ^ presents the gauge choice, the tensor-scalar equations 



derived in this gauge (Sec. 2.1) and the numerical method used for integrating them (Sec. 2.2), as well 
as the performed tests (Sec. |2.3D . Section ^ is dedicated to the simulations which were done; namely, the 
physical models used ( |3.1[ ), results of the simulations ( |3.2| ) and a discussion on the resulting gravitational 
waves and their interaction with the detectors (3.3). Finally, Section ^ summarizes the results and gives 
some concluding remarks. 



2. Equations and numerical issues 

As in Paper I, we have chosen the Radial Gauge and Polar Slicing (Schwarzschild-like coordinates), the 
metric being diagonal: 

ds"^ = -N^{r,t)dt^ + A'^{r,t)dr^ + r^{de^ + sin^ 9d(f>'^), (5) 

where N{r,t) is called the lapse function. The metric g*^ will often be described by the three auxiliary 
functions v{r^t), m{r,t) and (^{r^t) defined by 



iV(r,t) = exp(i.(r,i)). 
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A{r,t) - 1 ^ 



All coordinates are expressed in the Einstein frame, and asterisks are omitted. However, "physical" 
quantities will often be written with a tilde, that is in the Jordan-Fierz frame. The matter is supposed to 
be modeled by a perfect fluid: 

T^i- = (e + p)u^u^ + pg^^ , (6) 

where is the 4- velocity of the fluid, e is the total energy density (including rest mass) in the fluid frame 
and p is the pressure. 



2.1. Equations in RGPS gauge 



As far as the hydrodynamics are concerned, we will use a set of primitive variables {ns, e, U}, which 
are respectively the baryonic and the total energy densities in the fluid frame, and the fluid radial velocity: 

A u'' 



U 



r = (1 — U^) ^1"^ being the Lorentz factor of the fluid. We also use the variable: 

^--To" = r'(e+p)-p, 

and three "scalar-field" variables: 

Considering the conservation of the stress-energy tensor and of the baryonic number (Cf. Paper I), one 
can write the hydrodynamical equations for evolution in a hyperbolic formQ: 



du 19 

dt dr 



(7) 



where u = {Z)*,/i*,r*} is the vector of unknowns (evolved quantities), defined as 

D* = a'^{ip)ATnB, 
M* - aHv)[E + p^U, 
T* = a^{ip)E-D*. 



(8) 



^the equation for the density is not strictly in that form, see hereafter, at the end of the section 
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The associated vectors of fluxes f{u) and sources s{u) being (we define q-^ = 47rG'*/c*): 



fD' = UD* 
fr' = i-i*-UD*. 



(9) 



and 



= NAiy{^i*U -T* - D*)[2q^rp 
3a{ip)p*ri 



A 



2rfi*7]^ - -E (r* + £>* + p*) (1 + U^) 



(10) 



{t* + D* + p*) NAr {{1 + U^)i;r] + UE} -sd-- 



With this formulation, we easily recover the set of hydrodynamical equations in General Relativity, in the 
limit {if — > ipo, a{(f ) — > 1), as described in Romero et al. (1996)| . Here, however, the evolution equation for 
the relativistic density (D*) is not strictly in the form (^, but writes: 

dD* a{ip) d 



dt 



j,2 Qj, 



,N 
A- 



(11) 



The system of equations is completed by an equation of state p = pijiBi e), which will be detailed in 
Sec. BJ, and the field equations for the metric and the scalar field: 



dm 
dr 
dv 
dr 

ay 



—r^ (S + g.(r*+D*)), 
2 V47rr2 



dCdip\ dCdip 



(12) 
(13) 
(14) 



For more details on the derivations of this system, see Gourgoulhon (1991), Paper I and Romero et al 



(1996). The key point is the choice of the conserved quantities (|^) which, in tensor-scalar theories, keeps 
the system in a strictly hyperbolic form (|^). 



2.2. Numerical integration 



As it will be discussed in Sec. 3.1, shocks are expected to appear in the simulations. It is then clear 
that the numerical integration must be able to handle such discontinuities. The code used in Paper I is 
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based on pseudo-spectral techniques and so, is not well adapted for this kind of problem. However, in one 



dimension Bonazzola & Marck (1991) have adapted the spectral methods to the case where a shock was 
present. This technique will not be used here, since we are interested in generalizing (in a future work) 
techniques developed here to two- and three-dimensional calculations. All the details of the numerical 



techniques will not be described here, one can refer to Romero et al. (1996) and Marti (1997) for the 
Godunov-type (High-Resolution Shock-Capturing, hereafter HRSC) methods and to Paper I and references 
therein for the spectral methods. 

The idea is then to use two different numerical grids to describe the same physical space (using the 
same coordinates). These two grids allow for two different numerical techniques for solving the (^ - 
(p^), governing the evolution of the matter and fields. The HRSC methods, on a grid with a rather large 
number of points (typically a few hundreds), are used to solve the hydrodynamical hyperbolic system (j^), 
for they are the most efficient in handling discontinuous quantities, such as density or velocity after the 
bounce. Spectral methods are used to describe gravitational fields (g*^ and ip), which arc smooth, as well 
as their first derivative. In this work, for the equations of the Sec. |2.l| only the wave equation ( p^ ) has been 
integrated by those means, the two constraint equations, (|l^) and (|l^), giving the metric are very simple in 
the RGPS gauge and were integrated by usual finite-difference method. 

Last but not least, the numerical information of the various fields has to be passed from one grid to the 
other paying special attention to minimize the numerical noise, not to destabilize the integration schemes. 
Getting from the spectral grid (spectral representation of functions as a truncated series of Chebyshev 
polynomials) to the HRSC one shows no difficulties, since a function, represented by its Chebyshev 
coefficients, can be evaluated at any point of its definition interval with a high precision as a (truncated) 
series and using the properties of Chebyshev polynomials (see e.g. Bonazzola et al.(1997)| ). Nevertheless, 



the interpolation in the other way (from HRSC grid to the spectral one) has to be "smooth" enough in 
order not to bring high frequency terms. The method used here meets these requirements, since functions 



are interpolated in the smoothest way, introducing no spurious information (Bonazzola (1998)). The 



interpolated points are such that the resulting amplitude of the function /, calculated as 



' interval 

is minimized. This ensures that the Chebyshev coefficients of the interpolated function decrease fast enough 
and the overall series converges. 

The resulting code (RSSM in the next,RS stands for Riemann Solvers and SM for Spectral Methods) is 
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then just the combmation of an HRSC module, integratmg in time the hydrodynamical variables {ris, e, U} 
and resolving the constraint equations giving m and v. The source term of the wave equation is then 
evaluated on the HRSC grid and interpolated to the spectral grid. Then, the wave equation is integrated 
to get the scalar field, which is computed, on the HRSC grid from its Chebyshev coefficients. The typical 
number of points used for the HRSC grid is between 300 and 500; for the spectral grid it is between 65 and 
129. 



2.3. Tests of the code 



The HRSC part of the code comes directly from the code presented in Romero et al. (1996) , in which 
numerous tests are shown. The spectral part, solving the wave equation, has been presented and tested in 
Paper I. The interpolation procedure has been tested on various analytical functions. The combined code, 
as a whole, has been checked by setting the scalar field to its asymptotic value (po (as well as the coupling 
function a{(p) (Q)), which gave us the general-relativistic results. On the other hand, putting the pressure to 



zero enabled us to retrieve the dust collapse of previous studies ( Shibata et al. (1994) , ^checl et al. (1995 



and [Harada et al. (1997)| ) 



Finally, we tried to follow the collapse of a neutron star to a black hole, to compare with previous 
results from Paper I, but with this combined code. The initial model taken is a neutron star with a 
polytropic equation of state (see hereafter eq. |l^). Since HRSC methods do not achieve the same precision 
as spectral ones, the collapse has to be triggered in a artificial way: decreasing the adiabatic index 7 or the 
K coefficient of the equation of state: 

~/~ \ ~ ~ , z^^-O^lB f ^~bY 

e(nB) = nsmB + K — 

7 - 1 V "-0 / 

riB ^ 



p = KhorhB — ■ (15) 
KnoJ 

ifiB is the mass of one baryon (= 1.66 x 10~^^kg) and no — O.lfm^"^. The resulting scalar waveform (see 
Sec. p^ ), computed by both codes separately is displayed on Fig. |l|. The discrepancy is quite important at 
the beginning of the collapse, where the instability has to develop and the HRSC methods are less accurate. 
In addition, the boundary conditions imposed on the hydrodynamical system are not rigorously the same: 
as mentioned in Paper I, the spectral code uses a comoving grid, and there is only one boundary condition 
to impose which is chosen to be p-^?'^?7|^_^ — 0; the HRSC methods use an Eulerian grid and impose 
three boundary conditions, chosen to be /c* = Z^. = fr' = at the boundary of the grid. This can 
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Fig. 1. — Comparison between the pure spectral code and the combined (HRSC and spectral) for the collapse 
of a neutron star, of baryonic mass 2.2 Mq and with ao = 2.5 x 10~^, (3o = —4. Time origin is arbitrary. 
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influence the beginning of the coUapse, but after that, the evolution is fast enough and such discrepancies 
disappear. 

3. Simulations 

As it has been mentioned above, the physical scenario we are considering is that of the formation of a 
neutron star. 

3.1. The model 

One can think of two ways to come to a neutron star: the collapse of the iron core of a massive star at 



the end of its evolution, i.e. the standard mechanism of type II supernova (Shapiro & Teukolsy (1983)), or 



the accretion induced collapse of a white dwarf which, under some conditions, can be non-explosive ( |Cana] 



(1990) ). In any case, the degeneracy pressure of the electrons can no longer support the gravity of the core 
and the collapse is unavoidable. When the central density becomes comparable to that of nuclear matter, 
the strong interactions between nucleons stiffen the equation of state. A strong shock is then generated 



according to the well known model for prompt mechanism of supernova (see Miiller (1997) for details). 



Since we are interested in the emission of monopolar gravitational radiation in the framework of a 
tensor-scalar theory of gravity, we have simplified the micro-physics concerning the equation of state (see 
below) and the neutrino transport (which has been neglected). The equation of state we consider is that of 
a perfect fluid for which the adiabatic index 

dlnp 



7 = 



(16) 

adiab. 



d In fl B 

is not constant throughout the evolution of the collapsing star, but is a priori dependent on the density, 
temperature and chemical composition of matter. During the collapse, the adiabatic index is supposed to be 
lower than 4/3, but when the strong interactions between nucleons become dominant over other pressures, 
the adiabatic index grows rapidly above 4/3 and even 2. We can also consider that all the matter in the 



core follows roughly the same trajectory in the density-temperature space (see Arnett (1977) and Van Riper 



(1978) ). In a first approximation, the adiabatic index thus depends only on the density: 

7("s) = 7min + S (logio(p) - logio(pbounce)) (17) 

where p = rhsnB, and (pbouncc, S) are two parameters of the equation of state. There are other simplified 
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versions of the equation of state, which try to minimize the complex physics involved (see e.g. Miiller 



(1997)). In this work we use (|l^)and we have taken 7,„in — 1.33 in all our runs. Table 1 summarizes the 



different parameter used for the models studied in this work. Models A, C and E use "standard" values for 
the equation of state parameters, whereas models B and D will produce a "stiffer" bounce and allow to 
explore higher central densities. This second type of model, even if it seems less realistic, allows for more 
relativistic regimes (the resulting neutron star being more compact), with stronger shocks. The initial state 



is the same as in Romero et al. (1996): a white dwarf of 1.39Mq, which is the maximal mass for such 



objects, as derived in Ibanez (1984), Models A and B are the same as in Romero et al. (1996), and will be 
used as general-relativistic calibrations, to compare with tensor-scalar results. 



3.2. Results 



We briefly show the results of general-relativistic simulations (models A and B of Tab. 1). We retrieve 
the results of Romero et al. (1996)| , but with the combined code, which is evolving a constant (equal to ipo) 
scalar field (with a null coupling function (^)). Figs. ^ and|| gives the velocity profiles for both simulations: 
both show the characteristic "V" shape during the infall epoch (labels 1,2 and 3 for Fig. ||, and 1 and 2 
for Fig. H). When the density near the center becomes larger than ptounce (see Figs. ^ and ||), the infall 
is stopped and the matter coming from outer layers is stopped and a shock appears propagating outward 
(labels 4 to 9 of Fig. |2| and 3 to 9 of Fig. ||). Central density reaches its maximum when the bounce happens 
(Figs. 1^ and ||) and then settles to an equilibrium value. The evolution of the lapse is similar (see Figs. ^ 
and 0), and one can see that the final state of model B is more compact than A, even if the value of the 



central density is not compatible with these which are usually admitted for neutron stars (see e.g. Salgado 
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Table 1 
work. W 



Parameters of the equation of state and of the coupling function (m , for the runs presented in this 



; took 1 



1.66 X 10 ' kg • m 



200 



400 



radius [km] 

Fig. 2. — Snapshots of the fluid velocity U, as a function of the radius, at different moments of the collapse 
and bounce, for the model A of Tab. 1. Labels indicate time evolution (l:t=71.2 ms; 2: 76.3 ms; 3: 80.0 ms; 
4:81.3 ms; 5: 82.4 ms; 6: 83.5 ms; 7: 84.7 ms; 8: 85.8 ms; 9: 86.9 ms). 
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80 85 
time [ms] 

Fig. 3. — Evolution in time of the central density of the star, for the model A of Tab. 1 
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80 85 
time [ms] 

Fig. 4. — Evolution in time of the lapse at the center of the star, for the model A of Tab. 1 
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200 400 
radius [km] 

Fig. 5. — Snapshots of the fluid velocity U, as a function of the radius, at different moments of the collapse 
and bounce, for the model B of Tab. 1. Labels indicate time evolution(l:t=73.8 ms; 2: 79.6 ms; 3: 81.1 ms; 
4:81.8 ms; 5: 82.6 ms; 6: 83.4 ms; 7: 84.2 ms; 8: 85.1 ms; 9: 85.9 ms). 
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Fig. 6. — 



Evolution in time of the central density of the star, for the model B of Tab. 1 
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time [ms] 

Fig. 7. — Evolution in time of tlie lapse at the center of the star, for the model B of Tab. 1 
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et al. (1994)). 

We now turn to the tensor-scalar theory with the numerical evolution of models C, D and E of 
Tab. 1. Cases C and D have coupling function parameters that do not allow for the so-called "spontaneous 
scalarization" (see Novak (1998b)D , for j3o is not negative enough. Table ^summarizes some results of the 
simulations. The most important difference between these runs and the general relativistic ones comes 
from the gravitational radiation emitted during the collapse (see Fig. |ll| and |lj), which will be discussed 
hereafter. These two figures also show the evolution of the scalar field during the collapse, the differences 
between both being related to the fact that the case D is more relativistic than C. The resulting neutron 
star being more compact for D, the scalar field reaches higher values and its influence on the final state 
begins to be noticeable. In both cases, in the beginning (until t = 82 ms), the scalar field only "follows" 
the matter, producing the spike at t = 82 ms in both figures. Then, in the C case, the scalar field varies 
as the star's central density. Nevertheless, in the D case, the scalar field has not reached its equilibrium 
value yet, so it continues to grow. However, the differences between general-relativistic and weak scalar 
field runs (A and C or B and D) are very small, the only one noticeable is the difference in the kinetic 
energy of the matter reaching escape velocity, between runs B and D. This can be seen as if the scalar 
field deepens the gravitational potential (which will be confirmed by the study of the case E), by locally 



increasing the gravitational coupling function (see Damour fc Esposito-Farese (1992) ). The figures showing 
the evolution of central density, of the lapse and snapshots of velocity distribution are not displayed, for 
they are indistinguishable from those of General Relativity. 

The case E has parameters of the coupling function (^ allowing a priori for a spontaneous scalarization 
of the neutron star. The results shown in Tab. § suggest that has indeed happened: scalar charge and 
resulting gravitational radiation are (at least) one order of magnitude higher than for the C case. The 
figure ^ shows the fluid velocity profile after the bounce, which exhibits a series of oscillations behind the 
shock. These oscillations are no numerical noise since a change in the number of points in the grids gives 
the same result and, most of all, the frequency of these oscillations corresponds to that of a neutron star 



undergoing a transition to a "spontaneous scalarization" state ( Novak (1998b) ). It is particularly interesting 
to compare Fig. 4 in this reference and Fig. of this work. Let us recall that the numerical code used in 



Novak (1998b) is using only spectral methods and has been extensively tested; moreover, the particular 
simulation of transition to "spontaneous scalarization" has been verified by comparing the final data of the 
evolution (which represent a stable equilibrium configuration) with the results given by a static code. As 
it can be seen of Figs. ^ and |l^, the whole proto-neutron star oscillates, looses energy through monopolar 
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C\3 




200 400 



radius [km] 

Fig. 8. — Profiles of the velocity, as a function of the radius, for the model E of Tab. 1, before and after the 
bounce(from left to right and top to bottom). Times are the following (t=80.0 ms, 81.3 ms, 82.3 ms, 83.0 
ms and 84.6 ms). The oscillations behind the shock are due to the oscillation of the new-born neutron star. 
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Fig. 9. — Evolution in time of the density at the center of the star, for the model E of Tab. 1 
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Fig. 10. — Evolution in time of the lapse at the center of the star, for the model E of Tab. 1 
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43 


7 
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0.38 
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0.38 


0.64 


rdip [m] 








27 


100 


800 


Srad [1O-3M0C2] 








1.3 • 10-3 


3.2 • 10-2 


0.8 



Table 2: Various quantities measured during the numerical evolution of models described in Tab. 1. Mg 
is the gravitational mass of the resulting neutron star, to its scalar charge (defined as the monopolar part 



of the static scalar field, see Damour & Esposito-Farese (1993) for more details) and i?_3 is such that 
/5(i?_3) — 10-3p(r — 0). Skin and E^ad are the kinetic energy of the matter reaching the escape velocity 
when the shock comes to r = 600 km, and the radiated gravitational energy as a monopolar wave. Pmax is 
the maximal value reached by the central density A^min the minimal value of the lapse and rdcp is related to 
the amplitude of the scalar gravitational wave by (|l9|). 
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radiation and settles down to a strong scalar field state. We also have the confirmation that the scalar field 
deepens the gravitational potential well: comparing the run E to A and C, one sees that the neutron star is 
more compact and that the ejected matter has less kinetic energy. 



3.3. Gravitational waves 



We here recall the results of Damour &: Esposito-Farese (1992) and Wagoner fc Kalhgas(1997) 



concerning the interaction of a scalar wave with interferometric detectors. Far from any source of 
gravitational radiation and in a weak field, the metric writes: 

ffM- = a\^o) (^7?^. + [hl^ + 2aoFr^^,] + 0{l/R^)^ ; (18) 

where rj^u is the flat (Minkowskian) metric, h*^,Xt ^ R/c) and F{t — R/c) are the \/R components of the 
Einstein metric 17*^ and of the scalar field The Riemann tensor of the physical metric can be written (at 
the first order) as the sum of two terms: 

which is the spin-2 contribution to the wave {TT denotes the transverse-traceless part) and 



Rc^ \ dt^ dx^dxi ^ 
It is thus clear that monopolar waves interact with the detectors as well as quadrupolar ones. It is known 
that if one writes the geodesic equations in the transverse-traceless gauge, one sees that test masses have no 
coordinate change when a quadrupolar gravitational wave passes. But the photons measuring the distance 
between these two masses are affected by the wave. The opposite holds for monopolar waves: the photon 
travel is not affected (since they follow null geodesies, they cannot differentiate between Einstein and Fierz 
metric), but proof masses undergo coordinate change in the TT gauge as the wave passes. 

The right quantity to compare with the quadrupolar wave amplitude at a distance d from the source 
hq — hjj'^ /d appears to be: 

hs{t) = ^a^iipo)aoFit), (19) 



where d is expressed in meters. This is the reason why we have plotted the function F{t) in Figs. |ll[ - |l3| , as 
in Paper I. These three figures represent the three type of waveforms we have encountered. In (Fig. [Tl] ) the 
scalar field just "follows" the hydrodynamical evolution, having almost no effect on the hydrodynamical 
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time [ms] 

Fig. 11. — Resulting monopolar gravitational waveform (as a function of time), for the model C of Tab. 1 
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Fig. 12.— 



time [ms] 

Resulting monopolar gravitational waveform (as a function of time), for the model D of Tab. 1 
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Fig. 13. — Resulting monopolar gravitational waveform (as a function of time), for the model E of Tab. 1 



- 27 - 



evolution. In (Fig. ^2|), it acts a little on this evolution and undergoes raise even after the bounce occurred. 
Finally, in (Fig. |l^) the star undergoes a "spontaneous scalarization" and the waveform is very close to 
those shown in Novak (1998b) . The parameters on which the waveform depends are essentially the (3q 
(for the spontaneous scalarization to appear) and the final compactness of the star, which depends on the 
equation of state used ([l7|). The more compact the star, the stronger the scalar field. 

The results of Paper I, concerning the dependency of the amplitude on the scalar field parameter ao 
were retrieved; in the case of spontaneous scalarization the wave amplitude ( |l9| ) goes as ao, and as ciq 
otherwise. Since this is the amplitude which interacts with the detectors and since it contains at least one 
factor ao (coming from the interaction of the wave and the detector, in weak field regime), which is strongly 



constrained by solar-system experiments (Williams et al.(1996) and Damour & Esposito-Farese (1996)), 
the result may seem disappointing, if spontaneous scalarization is excluded (as suggest recent results from 



binary-pulsar timing, see Damour & Esposito-Farese (1998)) and if one compares to the amount of released 
energy (Cf. |[ However, if we consider the case of the model D: taking ao = 3 x 10^^, which is the maximal 
value compatible with the solar-system experiments, taking the detectability limit of VIRGO^ (Caron 



et al.(1997) ) 3 x 10 Hz at 800 Hz (which is the characteristic frequency of our events), one sees that 
a collapse occurring closer than 10 kpc from us (half of our Galaxy) could be detected. 



4. Conclusions 

We have successfully studied the gravitational collapse leading to a neutron star, within the framework 
of a tensor-scalar theory, and this is the first study of a gravitational collapse, with the presence of 
shocks, in tensor-scalar gravity. We have shown that, if the parameters of the theory are such that the 
non-perturbative phenomenon involving the scalar field (the so-called "spontaneous scalarization" ) can be 



excluded, as seem to indicate results of Damour & Esposito-Farese (1998), then the scalar field has a little 



influence on the hydrodynamical evolution of the collapse. If the case where spontaneous scalarization can 
occur, the run indeed indicated that the resulting neutron star had a very strong scalar charge and that it 
underwent strong characteristic oscillations. This gives us high confidence in our combined code. As far 
as gravitational waves are concerned, we have shown that, with ao and f3o constrained with solar-system 
tests and binary-pulsar timing, such a collapse could be detected up to 10 kpc from us. From the opposite 



■^at the considered frequency, the sensitivity of LIGO is comparable 
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point of view, a non-detection of any monopolar wave, when a supernova occurs (when we have neutrino 
and/or electromagnetic signal), could give constraints on the parameter space of the tensor-scalar theory. 
The computed waveforms appear to be quite different from those of Paper I, meaning that, if monopolar 
waves from gravitational collapses could be detected, one could in principle differentiate between each 
type of collapse: Type II Supernova or collapsing neutron star. As a byproduct, we have developed and 
tested a new type of code, combining spectral methods and Godunov-type techniques. We hope to proceed 
to more complicated problems (in more than one dimension) in General Relativity (that is, without the 
scalar field), using, at each time step, HRSC methods for solving numerically the equations governing the 
evolution of matter (the Hydro-part) and spectral methods for solving numerically the equations governing 
the gravitational field equations (Einstein-part) just for the distribution of matter given at that time step 
by the Hydro-part. The experience gained in present work seems to be stimulating for such new strategy in 
Numerical Relativity. 

We are indebted to Silvano Bonazzola for providing us with the "smooth" interpolation procedure, 
as well as for very useful discussions and reading of the manuscript. JN acknowledges financial support 
connected with a post-doctoral fellowship LAVOISIER from the French ministry of foreign affairs. This 
work has been partially supported by the Spanish DGES (grant number PB97-1432). 
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